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АВ5ТКАСТ 


A study is made of the circumstances under which there will 

be incipient buckling in an infinite plate containing a circular 
hole and subjected to uniform uniaxial tension at a great distance 
from the hole. This buckling is a consequence of the fact that 
near the hole and in line with the direction of the applied tensile 
loading, there are areas in which both principal membrane stresses 
are compressive. Values of critical tensike stress are given, 
together with expressions for the corresponding lateral deflection 
patterns. The results are in the form 

S = ee у. 

cr 
where E is Young's modulus of elasticity, t is plate thickness, a 
is hole radius, and K is a buckling coefficient having the following 
values for the first four modes: 7.65, 7.80, 11.7, 12.3. A 
comparison is made with some experimental determinations of Danis 


for plates of finite width. 
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1. INTRODUCTION 

The problem of determining the stress situation near a hole 
in a plate under tension (see Figure 1) has a well known solution 
which is of great practical importance. 

The original work on this problem was done by Kirsch (19), 
who solved the problem of the stress distribution in an infinite 
plate pierced with a round hole and subject to uniaxial tension 
at a great distance from the hole. Howland (17) treated a similar 
problem with an infinite strip of finite width. Savin (5) and 
Hogan (15) have compiled solutions produced by other authors for 
variations in hole geometry, as square, triangular, etc. for both 
finite and infinite plates.  Lur'e (3) has solved the problem with 
a hole in a shell of uniform curvature. 

An important result obtained from these solutions is that the 
stresses near the hole become much greater than those imposed upon 
the plate at a distance from the hole. This has given rise to the 
wide publication of many stress concentration factors which are 
important in design work and which relate the maximum stress to the 
nominal stress at a distance from the hole. 

Another feature of the stress situation, which is not so widely 
recognized, is that certain areas of the plate are under uniaxial or 
biaxial Popes sone although tension is imposed externally. These 


areas of compression, giventhe proper conditions of plate thickness 


* These terms will be employed herein to mean that one or both 

principal stresses are compressive. Specifically, "uniaxial compression’ 
means that one principal stress is compressive and one is tensile; 
"biaxial compression" means that both are compressive. 
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and external tensile load, will buckle, possibly causing failure of 
the plate. Houghton and Rothwell (16) have shown that in plates 

with holes having a very small gross curvature the stress concen- 
tration factor can increase appreciably. It is reasonable to assume 
then, that even though local distortions of a plate may be acceptable, 
such distortions could increase the stress concentration beyond that 
expected for a flat plate. Failure of the plate by yielding could 
then occur. 

There is no known investigation of this buckling problem, other 
than an experimental investigation by Danis (14), to determine the 
critical stress for the buckling of a finite plate with a round hole. 

This paper presents a theoretical investigation of the incipient 
buckling problem for an infinite plate containing a round hole. This 
problem is of interest in itself and for the practical reasons given 


in the preceding discussion. Results are presented on page 33. 
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2. PROBLEM DESCRIPTION 


Consider a plate infinitely extended in two directions and 


pierced with a hole, as shown in Figure 1. 


S 





Figure 1 

x and y are Cartesian coordinates, r and 9, cylindrical coordinates. 
The z axis, not shown, is perpendicular to the paper at the center 
of the hole, and extends in the positive direction toward the reader. 
Consider all quantities to be in units of hole radius, a, which 
essentially reduces them to dimensionless quantities, with hole radius 
equal to one. 

The plate is placed in uniform tension, S psi, at infinity in 
the x direction. The solution of the stress distribution problem for 
this plate, as originally solved by Kirsch (19), and clearly presented 


by Timoshenko (7) is: 
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This set of equations is henceforth referred to as the Kirsch 
solution. 

It is observed that the maximum tensile stress is 3 times the 
stress applied at infinity, and occurs at т = 1, апа Ө = +T/2. 

Ihe maximum compressive stress, equal in magnitude to the applied 
tensile stress, occurs at r = l, and 0 = +7Í, From this it can 
be seen that if buckling were to occur in the plate, it would be 
most apparent near the hole along the x axis. 

А better understanding of the buckling problem is obtained by 
considering the minimum principal stresses (most negative, largest 
compressive). From the Kirsch solution, the minimum principal 
stresses are found and shown in Figures 2 and 3. 

Figure 2 shows one quarter of the plate, with the hole centered 
at the top left. The stress distribution is symmetrical with respect 
to the x and y axes, This picture was prepared with consideration 
only to whether the stress situation at a particular point is biaxial 
tensile, biaxial compressive, or uniaxial compressive. 

Figure 3 shows the same section of plate as Figure 2. This 
presentation maps the relative values of the algebraic minimum 


principal stress at each point. 
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ч, чч, “ч, че, ч, ч, ч, ч, “чч ам, ч, ч, ч, б, ч, ч, ы, ч, “ы, бы, ы, ы, бы бы, бы бы ы ы ы ы а 
a e ы “ч че, Чч, Чы, “ч ч, ч, ч, “ч, чы, ч, ч, ч, ч, е, б, ы, ы, “ы, м, ы, ы, ы, ы, ы “ы ы ы л 
НУ 
Но 
BE ч “ч ч, 6 
CMM, А N N N =, “=, чу Зе, Эч, “ч “эъ, ч, е, ч, ч, б, ы, б, ы, ы, ы, а, ы, бы, ы, ы, ы, 
МУ ч, “ч, чч, Че, чы, ч, ч, ч, ч, “ч, а, б, ы, N N N N N бы, ы, бы, Мы, бы бы, “ы, ы, ы, бы ы, э 
w... “ч, “ы, “ы, ы, ч, ы, ч, ч, ч, а, а, Чы, бы, Мы, б, ы, ы, ы, ы, ы, ы ы, ы, ® 
мы ч “ч “ч Эч, “ч “ч, e, ч, ч, “ч, ы, ч, ы, ч, ы, ы, чы, ч, ы, б, а, ы бы, чы ы ы бы, ы ы 
N ч Чч, Чч, “ч, Же, ы, “е, ч, ч, ч, а, ч, ы, ы, а, ы, а, б, ы, бе, а, ы, ы, ы ы ы Мы, ы бы 
Мы ч, Чч, ч, ч, Же, ч, ч, ч, “ч, ы, б, ч, ы, ы, ы, ы, ы, б, ы, ы, ы, бы, ы, ы, чы ы Мы, ы ы ө 
ч, “чч, чч, “ч, ч, ч, ы, “ч, ы, ч, ч, BWA ы, ч, а, ч, ы, ы, ы, б, ы, ы ы, а, бы ы, чы б ы, ы, ы ы, ы ы ы ы быз 
SO “ч “ч “ч, “ч, Че, “ч, Мы, ч, ы, ы, ч, ы, ы, ы, ы, ы, бы, ы, ы, ы, ыы, аы, аы, ы, бы N аы, Чы Чы, ыы, бы, ы ы ы ы чы 
ыч “ч, “ч, чч, ч, ч, ч, “ч, е, ы, Ме, ч, ч, б, ы, б, Жы, а, ы, а, ы, ы, ы, ыы, бы, а, ы, бы был ы, быы, ы, ы ы, ы ы ы ы ЭЪ 
K w... N... ҸҸ аы, б, ы, ы, ыы, бы, аы, ы, Жы, бы, ы ы ы Чы ылы 
SNS N N N N N N N N N N N N ч, ы, ч, чы, "ы, ы, ы, ы, ы, ы, м, Мы, а, ы, ы, ы, ы, ыы, ы, Мы, ы, ы ы, Чы, ы, ы, 


TENSION 
FIELD 


N SS a ч, б, ы, а, ы, ч, ч, а, зы, аы, ы, ы, аы, бы, ы, FR, Чы FRU, бы, ы ы, Чы ы ы a ы, Чы, ы ® 

Мы ч, “= ч, “ч, Чч, ч, “ч, “ч, ч, ч, ы, ы, ы, ы, ч, ы, ы, ы, м, бы, ы, ы, бы, ы, ы, ы ыл, ы бы бы PR, RH RR RR 

TAS, RC, “ч, ы, ч, ч, ч, ч, ы, ы, ч. а, ы, бы, ы, ба, ы, ы, аы, бы, ы бы, бы, ы, ы, ы, ы бы б ы, Жы бы, ы ы ы N 

N ч М чч ч, чч, ч, Че, ч, ч, "ч, ч, ч, “ч, чы N аы, чы, "аы N N NN N аы, ыы ы, бы, Мы, быы ы ы бы ы Чы ы ы бы ы ы Мы ® 

N SN ON NON N N N N N N ON N N N N N N N N ы, ы, аы, М, бы, бы, Мы, бы аы, Мы бы, бы, ы ы ы ы ы ы ы ы ыы 9 

ММ ч “ч чы, ч “ч, ч, че, “е, ч, Че, ч, ч, “ч, ч, о, ч, чы, чы, ы, ы, ба, а, ы, ы, Мыл ы, ы, ы, бы, МЫ Мы ы, ы, Жы, бы, ы ы ы ы a, + 

ч А е, “е, ч, “е, N N е, N N N ы, ы, ч, ч, ы, ы, ы, ы Аш ы, аы, ы, ы, ы, аы, ы бы, ы бы бы ы бы бы, Жы ы ы Мы HR RR 

a “ч, “ч, “ч, UM, ч, ы, ч, ч, чы, ы, б, ч, б, а, ы, чы, ы, чы, ы, бы, ы, ы, ы, бы, ы ы ы, ы ы ы, ы ы, бы бы, Чы Чы бы ы чы © 
TQ, S, a, Чч, Чч, ч, e, SN NN N ON N m, ы, а, а, чы, ы, ы ыы, ы, бы бы, ы ы ы ы ы ы ы бы бы ы чы бы ы ы ы ы ы 
ККК 
NSS N ККК КМ ы, ы, ы, ы, ы, ы, ы ы чы ы ы бы ы ы ы ты ы ы чы ылы ® 
МОУ ч, =, ч Чы, Че, ч, ч, ы, ы, е, ы, е, ч, ы, М, ы, а, ба, б, Жы, ы, а, бы, ы, Мы, Мы, ы ы ы ы ы Мы, ы ы ы ы Чы ы ы Мы бы ° 
N NN N N чч, е, ч, e, ч Жы, ч, “е, Же, ч, Же, “е, е, ы, ы, ы, %, ы, ы, ы, ы, б, ы. бы, ы ы Жы Мы, ы Чы ы бы ы бы ы Мылы 
o “ч ч ч, “ч, ы, ч, ы, б, ч, ы, е, ч, ы, ы, чы, ы, ы, ы, ы, бы, ы, ы ы, ы ы чы ы ы ы ы Мы ы ы а ыы PR М 
ММ MS, =, Чч, “ы, Чы, ч, е, ы, ч, ы, ч, ы, ч, ы, Bea аы, ба, Мы ы, ыы ы, ы, Жы бы TR, А ы ы ы чы бы ты ы ы ы чы бы ты Мы Чы ы лы ылы бы ы 
~~. SS SSAA AN ч, ы, ы, е, Мы, ы, ч, ч, Ма, а, б, б, ы, ы, ы, ы, ы, ы, КМ 
AS. SSA RA RAR ы, ч, ы, ч, е, ы, ы, ч, ы, ч, а, Чы, c, T RM, а, ы, Мы, ы бы, ы Мы, ы ы ы бы бы ы чы бы ы бы бы Чы ы бы Чы Мы ыы бы ы 
ыч чч UM, Чч ч ч, е, ч, ч, “ч, ч, Эъ, ч, ч, б, чы, ы, а, бы, ы, а, ы, e, аы, аы, б, ы, Жы б, чы бы ы ы ы Чы ы ы чы бы лы ы ылы ы ы быы ы ы Чы лы Мы Мы 
NN ЛА Чч Cc, ч, ч, ч, Чч, е, “ч, Мы, ч, е, чы, ы, ы, ы, ч, ы, Чы, ба, а, аы, б, ы ы бы, ы Мы ы ы ы А ы Зы ы ЧЫ О ЫЫ ЧА лы ы ы чы быы Жы ы ы Мы 
N ы “ч, чы, “ч, чы ч, ч, ч, ы, ы, ч, ы, ч, б, ы, ы, ы, аы, Мы, а, ы, бы, аы, Мы ы, бы, бы М 
ММ ч “ч, “ч, е, ч ы, ч е, “ч чы, “е, ч, ы, бы, ы, Чы, е, ы, б, ы, аы, ы, ы, ы ы, ы, ы, ыы, ы, ы, ЖЫ ККМ чы. 
NS N N N N N N “ч, ч, ч, ч, ч, Же, е, б, а, ы, б, ы, а, ы бы, Жы, а, бы, ы, ы Мы ы, ы, ыы ы ы А Чы бы бы ы ы ы ы a ы Жы Чы ы Жы Чы Мы ылы ты ы 


X — DENOTES BIAXIAL COMPRESSION 
/ — DENOTES UNIAXIAL COMPRESSION 


МОМ е ч S, f, ч, ч o a a “ч, ч, ч, “ч, ч, y, “ы, o, ы, ч, ы, ы, Эъ, а, a ы, б, ы, ы, Жы аы, ы ы ы, ы аы, a a ы, аы, a a ы, ы, Ж ы ы Чы ы, бы ы ы бы ы ы Ча ы, Чы, a ы, ы, 
.S SSS 55555 55.5 555 858 8 S55 8 4 84 S58 4 4 55 5 5858 5 8 8 8 85 8 8 858 858 8 8858 58455 85555 55 5455 
SA SSS 58 5554 8 5 8555 8 2585 4585 S545 85 5858 4 55 85 8. 8 8 8 858 58 28 858 45 55855555555, 88585555 
Н a, 


XQ, “ч ч “ч чы, Ж, ч, ч, ы, “ч ч, ч ч, чы, ы, RS, TR, f, чы, ы, TR, ы, Ж а, ы, аы, аы, аы, ы, ы, ы, ы, ыы, аы, аы аы, sq аы, RN, ыы Жы, ы, ы, аы, ы, ы, бы ы, ы, ы, бы, ы, Чы, бы, ыы, ы, а, ы, o 
w... T... N N N N N N N N N n N n а, ыы ы, ы, ы, ы аы, ыы ы, чы ы, ыы, аы, ы ы аы, Чы Зы, ы ыы Жы ыы, аы, ыы, аы, аы ы ы, бы Чы бы, ы, ы ы, ы ы, ы, ы ыы, ы, ыы, Чам, a, 
w... чч, ч ^ч, 58555 55 55.5 284 2 5 5 4 5 ы, ы, ы, ы, ы ы, ыл 8 5 8 аы, ы, 8 4 58 84 8 8 2 84 4 8 88 858 5858 85 55 585 545° 


~~. SS SS S54 S58 5564, 25 8 8. 25 55. 25 25 85. S58 2 58 588 85 4 8 5 55 858 858 858 8 58 58 88 5 85 5 558 5 5 55 854,555 


SSSA 554 554 S45 555564 55 5 855 5 5 55 558 5 555 55 58 55 58 58558 85 558.568 555585545 5545585555455 u — n 


QM, SAS 4S 4548 54 8 5 558 25 8458 4.58 8 85 85 4 458 4 8 85 85 85 858 85 45 8 588 4588584455555 5548 
~S SSA SS SSA SRSA ARR RS SRS RSS ASRS RS SRA SSS SSS 
TX M, M, S484 S858 ы, ч, 8 4) 58 8 8 4 8 8 85 858. 8 8 8 8. 858 85 4. 885 8855 88 ы ы, ы SASS 
SA. SARS SARS DAS WSEAS SSA RRR RS SS SRS SRS RSS SS 
AS А А Чч ч, SAAR AAS SSSR RSS SSS SSS SSS SSS SSS SSS SSS ы, Мы, ы 
~~. SAS SSS ASS SSS SSS SSSA RRS SSS SE SSS RASS SSS SSS a, 
SSS S54 5 4.88 8 5 8 4. 84 85 48 85 84544 8 585 555 55 4 8 58 85 554 858 8 85 5 58 85 88 885.4 55 85 5 854555 
+. SSS 45585 558 5 22 25 S45 25 5 855.558 S585 55 85 58 555 5 85 8 8 8 858 8 45884 8 8585 5858 55 558585555 
a O 458 558 8 854 58 55 45 402 2 5 8 858 2585 85 5 5 85 558 48 8 8 58 8 8 8 88 8 558 88 45 85 85 485 84554 5855 
SSS S548 S584 5.5 85 858 4 8 55 85 4 8 545 25. 85 85 558 58 5 55 8 8 85858 58 8 8 8 8 458 554 8585 858585555 
+A. SSA SRSA SSS SRS SSS RRS SSS ASS SS SS аы, ыы ы ыл ы ы ЧЫ ы Чы бы бы, бы, бы, Чы ы, бы бы, a ы, бы ы, ы, ы, 
ч ч чч NA A R n n n N N n N N n N 8 884 8 8 8 8. 85 84 88 8 88 88 8 4 84 558 ы бы ы ы 85545 
+S. SSR SOS SR NS RSS SS SSS SRS SS SS SSS SSS SS a a, 
+S SSS 5555 5454 55 5555 5 8 55 558 5 55, 85 55 85, 55 5585 5 85 858 85 88 55 858 55 55545 5545 5 455% 
ss... SSSA AS SARA SRS SR RSS RS RRS RR RS SSS SSS 
SMS, NN N a N 48 48 5 85 8 8. 4 8 85 8 858 8 588 48 8 258 8 ы 88 8 ы ы 4 85 86 85 8 5584585455 445 
SASS SS S555 5 55 85 55 5 585 2 5 55 258 858 S58 55 25 5 5858 2 85 28 85 858 858558455 5. 585 8585 8585 555 
SSS SSSA SAWSR RW SSR SW AVS BQ — — ` — ` — — — — — ` — — — ` — ыы, быы ы, ыы — — — — — — — — — ы, ы — — — — 
SSS S55 45 45 5 55 55 55 85 55 S58 85 8 5 85 8 5 88 85 84 85 254. 88 858 S58 8585 84 585 85858 58 8 8 585 5585 
Эч, “ч, “ч, ч ч ы аы бы, бы ыы, ы ыы ы N N Чы Мы аы ы аы Мы аы, Чы ы ы Чы ы быы ы Ма ы Чы бы ы ы ы Мы ы аы, бы, бы, ы ы Чы, бы, Чы бы, быы ы быы, быы Мы бы ы ы бы HN 
МО ч, ч, “ч “ч, ы, a ч, “ы, ч, ы, “чы, аы, a “ч, ы, “ы, ы, ы, ы, бы ы, ы а, ы, аы бы бы ы "аы ы ы, ы ы, ы ыы ы ы ы ы, ы ы ыы ы, аы, бы Мы бы, быы, бы, бы ы, бы, бы ы, ы, быы, бы, ы, ы, ы, 
ы “ч ч е “ч ч, чы, N N N N ч, ч, “чы, Эъ, ч, а, ы, “чы, а, ы, а, ы, ы бы, ы, ы, ы ы, ы, ы, ы бы Чы ы ы ы ы Чы ы ы ы бы ы Мы ы чы ыы ы ы ы бы, ы бы ы ы ы бы, ы ы ы ы ы Мы, 
К 
М 85 
К 
MSS, чч “ч, “ч, ч, чч, “ч, ы, a “ч, ч, ч, 8 5 58 54 8 5 25 2 8 558 8 8 55. 8 2 8 55 8 аы, ы ы ы ы ы 258 88558 58 588 858 45858 84864 85 85555 
a, 5 S54 45 854 5 5 55.55 55858 {84 88 8 258 855 54 858 845 5454 585 5858 8 28 55568 585884 555 585854 55 5 58584545 
SS... S58 55.554 55 85 8 858 8 чы, 8. 84 854 24 8 ы, ы, ы, ы, 8. 54 55 5 85845 8588 85 8 5554 858585554 05855 55 854545 
МОМ А “ч лы, “ч ө, “ч ы, “ч RR SAR RRS NS SRS HS ESSE SSS SSS 
N SN чч, ы, 55 845 858 8 5858 ы, 84 85 85 48 S58 5 855 5555 {8 4 28 585455 2584 5 85 25588 44 45464 58585854555 
МУМ М М чч “ы, ч М, “ы, “ы, Эъ, ы, ы, ы, “аы, N N N аы, ы, ыы, бы ы ы "аы ы, Чы ы аы ы, ы аы, ы, N ы бы ы, ы, ы, ы Чыл Чы ыы аы, аы ы, ыы, senses 
Sa ч чч ч, ч, “ч “е, ч, ы, ч, ы, ч, а, чы, а, ы, ы, "аы, а 864 54, 5 5858 855 54 86855 8 455 55 855 855555555 5 5555 5585845555, 
Мы ч, “ч, ч Мы, ч, ч, ы, чы, ч, чч, ы, ы, ч, ы, ы, ы, ы ы аы, “ы a ы, "аы ы a ы, nq Bw аы a {8 ы ыы, ыы ы быы, ыы ы, ы ы ы Чым ы, ы, чы ы ы ы ы ы ы ы, ы ы ы ы бы ы ANNAN 
NS чч ч чч, ч, Же, ч, ч, Эъ, y a ы, ы, ы, аы, ы, ы, чы аы ы, ы, ы ы, ы ы бы ы ы, a a бы чы N ы, Жы, ы a бы ы, ы, Чы быы ы аы ы N Мы, ы, a ы ы ы ы ы ы ы ы, ы, бы бы ы, ы ы М, 
Мы “ч “= Чч, “ч Ме, чы, М, Эч, “ч, “чы, ы, “ч, ч, ы, ы, а, чы аы, ы, ы, аы, аы, ы Мал, ы ы ы ы, ы ы ы быы ы 9, FR, FR, ы, 58 858585 45 858 588 455 84 285 8558555 
М ы ч, O Чы, “чы ы, “ч, ч, ы, ы, а, аы, ы, а, “а, ы, а, бы, ам, чы аы аы аы ы ы "ы бы Жы аы, ы, ы ы бы Чан ы n ы ы ы Мы Чы ы ы ы А ыы бы ы ы ы бы бы ы бы бы ы ы a HR 
SSS S555 6 55858 S58 5658 88 2 258 8 558 25 5585 2.8 S58. 25 25 258 258 25 85 558 858 864 85 5855 58548 585555 
o ч ч чч, ч Чч, ч, Чы, ч, ы, ы, чы, ы, ы, ы, ы ы, ы N N ы, ы, ы, чы бы ы ы ы ы ы ы, а ы ы, ы бы чы а ы, ы, ы ы ы ы ЧЫ Bw wna ws 6a sana 
~SS QS S525 55 64 25 85 S58 82 88 85 858 458 85858 84 5 5 55 S585 5 88 85 85 8 58 8 85 85 480.85 454 555 
SSS SSH... Мы, ы, ч, б, “ч, аы, ы, ы "ы аы, ы, Aq neq eqn ы Мы, а, ы ы чы ы ы ы Жы N ы Жы ы бы ы ы бы ы бы Чы бы ы ы Мы ыы 
ММ ч S55 55555. 45 8 24 S58 ы, ы 8 2 8 8558 85 85 8 5 ыы 5 5 58 8 558. 55 8 8 5 5 8 8 8 58 8 8 8 85 8 48 8558 584 8555, 
SSS 6 S64 85 558 S62 585858 2 S58 аы, ы, 2 4 2584 84 S858 868 858 S58 58 585 85 8 8.8 8 8848 64 858 884 8585 85585 85 55 55 
МОМ ч “ч ч чы, ч, ч е, BWW sq nq gg Qs Qe Bens Shes ыл ы ы ЧЫ ы ы ЖЫ МА Мы ӨЫ ы ы ы ы Мы ы ы Жы Чы бы ы ы бы ы ы ы ы ы бы ы ы. 
МОМ А ы чч Чч, Чы, ч, е, “ы, ч, ы, ч, ч, ы, ы, а, ы, ы, аы аы, а, "аы, а, ы, а, TR, ыы аы, ы, ы ы ы ы ы ы Жы ы ы Жы ы ы Мы ы ы ы ы ч ы МЫ ы бы ы ы бы бы бы ы ы бы ы ы ы Мы 
N SN SN N Чч ч, Зы, ч Чы, ч “чч, Мы, Мы, ч, чы, чы, ы, ы, ы, ы, ы, Ча, аы, ыы ы, а а, аы ы, ы, ы Жыл Чы ы ы, бы, Жыл бы ы ы ы ы Жы Мы ы ы ы ы ы ы ы, ы ы ы ы ы, ы ы бы, ы ы ы ы ы бы ы. 
SSS. SSSSL 55585856 25 252 ы, 2 58 858 84 58 55458 28 858 85 258 5. 558 45 5 85 88 8 8 84 2 8 54 4 88585 55 85455 
МЫМ ч “ч чч, чч, ч, ч, ч, ы, чч, Мы, VQ Qs. es. ыы, ы ы, аы, ы ы, ы аы ы ы ы ы ы ы ч ы бы ы ы na бы ы, ы ы ы ы ЧЫ ы Чы ы ы ы ы бы ы Чы бы, бы ы ы бы бы ы ы ы 
М Se чч чч, “ч Ме, Чы, Чы, “ч, ч, “ч, чы, аы, чы, чы, ы, аы, ы, а, ы, ы, Чы, ы, ы 258 858 5 858 85 858 8 258 858 858 58 88 58 858 4 858 885 5 45 4548 
SONS SS N N N N N N N N N N N N 22 2. 82 88 4 8 8 2 85858 8 5582 858 8 58 58 8 5 864 88 28 8 4 88 85885 864 бы бы, бы, ы, 85 80 4508 
~SS.S SSS S464 2585585658 858 2 8 8 28 28 28 8 554 8585 6.8 858 88 48 858 48 5 858 8 8 8 858 864 88 558 48 85855 55, 
.S SS 55858 8 56 8558 85858 8 85. “а, ы, 8 2 858 8 28. 258 858 48 5 5 8 85 258 854 85 28 8 8 8858 8585 4 5585 895 0 585 
МОМ ы SSS S554 84 55 58 5885 88 8 252 28 25 25 258 48 5 85 8 4 2 888 85 84 85858 85554 4858585 8 558 

МОМ Мы ч, “ч, “ч Чч, ъч, “ч, “ч, "ч, а, КМ 

w... чч, ч Чч, “ч, Чч, чы 2554 8 аы, 858 24. 5 8 44 2 88 8 858 284 8 8 85 85 8 88 48 4 8 8 8 88 88 84 58055 

~SS S855 5 S54 55 S58 85 88 22 8 4 2888 8 2 882 4. 28 8 4 545 84 SR 88 88 854 ы ы Мы Мы ы ы ы Өм 

aS SSS 55 5555 55 58 5858 S88 a 828 2. 85 8 88 858 ы 8 4 25 88 8 8 858 84 85 858 858 554558 

o o, qa Van nq nx as naan Мы ы ы ы ы 

SSS 55 S55 5 55 85 8.55 85 858 25 25858 45585 88 8555 55 5558 5858 55585 

SSS SS SSS SSS ы. ы, аы ы ы, ы, ыы, а М Мы ы ы ы Чы МЫ ы ы ы Жы, ыы ы ы ы ы ы ыы, 

SN SS e N n N N 5555 558 558 8 8 5 8 8 6.8 258 852 885858 8 8 858588 5 858 855858 

SSS 555 5. 585 5 5858 8 8 858 8 8 8 88 2 88 2.5 85 8585 88 85 8558 

Te, SS {55 55855 S58 554 858 5 4558 8 55 858 8 58 4 8555858 

~S S54. 4 4558 S555 58 5 8 858 4.5 84 8 55 585 55 58 5885555 

SSS SS S55 S55 55558 5 5558545 58 585 5855 5 55858 555 

a O O е, ч, ч, е, ы, ч, ы, ы, ы, аы, Мы, ы, ы, ы, ы, ыы, ы, ы, ы ы, Мы, ы чы 

Ме ч, “ч “ч, “ч Мы, “ч y a ч, ы, е, ы, ы, аы ы, ы, ы, ч N ыы, аы ы, ы, ы, Мы, 

SSS S556 545 458 4. 5858585858 5 5558 
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It can be seen from these figures that there is an area having 
biaxial compression. Accordingly, buckling is expected to take place 
for some sufficiently great value cf applied load. 

Numerous methods have been developed for solutions of plate 
buckling problems. In general, these methods involve either the 
satisfaction of the equilibrius equations or minimizing the strain 
energy of the plate (Energy Methods). A compilation of many solutions 
to plate buckling problems was made by Gerard and Becker (11). This 
is an excellent reference for beth method and practical solutions. 

One of the first energy methods was developed by Ritz (20), and 
independently, by Lord Rayleigh (4). This method is presented in a 
rigorous manner by Sokolnikoff (6). Budiansky and Hu (12) give a 
detailed example of the application of this method to a plate buckling 
problem. 

The Rayleigh-Ritz method involves determining the total strain 
energy in the plate as a function of an assumed buckled configuration. 
Equilibrium, and the stable buckled derlection, W(x,y), exists when 


this strain energy, U,, is minimized with respect to W(x,y). 


E 
fo do this, first, a family of functions, Фф , are assumed, which 


are dependent upon parameters A, , such that the buckled deflection, 


ke 
W(x,y), can be represented by: 
Ww(x,yY) = Plx,y, Ap Az) > 29 = 
where ф is such that the boundzry conditions on the plate are 
satisfied for any values of the А. 
The functions $ are substituted into the strain energy expression 
U, so that 


Jt = [fien | esa AA. EE ) dx dy (5) 


where U. is the integrated total energy in the plate. 
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In order to minimize the energy function, the following set 


of equations iS formed QU. 
کا‎ el О 
дА, 
gu 
о е (6) 
5А, 
If this set of equations has solutions E the function 


W(x 4) = P(x, y, A, Aaj is ) 
represents the buckled equilibrium condition, the boundary 


conditions of the problem being satisfied by $ for any values of the 


A... 
i 


2 


3. BOUNDARY CONDITIONS 

For the problem of an infinite plate with a hole, two boundaries 
are of concern. At infinity, deflection, slope, and curvature are 
zero. Around the hole, a free edge condition exists. 

The Kirsch solution is used, which already implies zero tractions 
around the hole, and uniaxial tension at infinity, so these conditions 
have already been satisfied. 

For a free edge to exist around the hole means that there is no 
shear, twist, or moment applied at the edge of the hole. Following 
Timoshenko (8), and considering the cylindrical coordinate system 


established in Figure 1, the conditions for a free edge are 


EA — (shear = 0) (7) 
v 

(М = = (twist = 0) (8) 

ce = O (moment = 0) (9) 


Kirchoff (18) has shown that the twist condition can be reduced 
to a statically equivalent shear force. The original requirements 
of zero shear and zero twist can be reduced to one requirement such 
that the sum of the shear forces and the static equivalent of the 


twisting moments equals zero. The conditions then become 


OMet | 


90 r Y= 


= О) (moment = 0) (11) 
r=) 


od) (Kirchoff shear = 0) (10) 
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My 








where Q, M, and M, are expressed as 
і тї 
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Mec? 20-94: AM - LA (14) 


When the appropriate differentiations and substitions are made, 


the boundary condition equations become 








(Kırchott 
ow ow aW | (2-39), 9W (3-2) zw o 
Кс a EIT ON E e ° -— = . == 
Ё "rar wor T vz “Sreez "үз КЕ gn ien. 
aw, 3[9wW, 1, 9?W|| _ Moment = 0) (16) 
gu. э | Әү +. кеш. 
r ү =00 


As previously noted, each member of the family of functions 
must satisfy the boundary conditions in the Rayleigh-Ritz method. 


For a choice of $ let 





202 А мм COS MO 
N=l M=0 Bum SIN MO 


where AM and Bm are arbitrary coefficients that are varied to 
minimize strain energy in the plate. 

The boundary conditions of eq. 17 are satisfied by this series 
since at r =&,W and all of its derivatives are zero. Conditions of 
equations 15 and 16 are satisfied by substituting equation 18 into 
these two boundary condition equations (see Appendix A). Such a 
substitution results in two equations, 15 and 16, with two unknowns, 
XM AND Y uM? for each conbination of N and M. A unique value of 


x 
X M and YM is found which allows each term of equation 18 to 


satisfy the free edge boundary conditions for any N and M chosen. 


From Figures 1 and 3, it is seen that there is symmetry with 
respect to both the x and y axes. This means that when W is written 


as a Fourier series in ® (eq. 18), the series contains either cosine 


* A term in this sense means the full expression 


Ana cos МӨ 
eee ient 





yt Bum SIN MO 
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terms or sine terms, but not both. Additionally, the requirement 
of symmetry around both the x and y axes means that only odd or 
even terms in M need be taken. Thus the Fourier portion of 


equation 18 reduces to: 


ANM cos MO. with M = 1,3,5,... 
or 

ANM соз МӨ, with M = 0,2,4,... 
or 

Bm sin MO, witheM-—-1,3,5,... 
or 

B sin MO, with M = 2,4,6,... 


NM 


This gives four separate expressions for W, one for each 
possible condition of symmetry. The full Fourier series could be 
taken, but if this is done, three quarters of the coefficients, AM 
or Bam are zero and considerable computer storage space is wasted. 


It is shown later than each of these series represents a different 


buckling mode. 
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4, ENERGY IN THE PLATE 

To this point, a family of functions W(r,9) (equation 18) have 
been developed, which satisfy the boundary conditions independently 
of the coefficients, Акм or ВАМ? and are adequate to represent 
possible shapes of the plate. Attention is now turned upon formu- 
lating an expression for the total strain energy, and then minimizing 
this expression to obtain the critical buckling stress. 

The energy in the plate can be divided into two forms. First, 
there is the membrane energy which is the energy stored because the 
mid-plane is stretched. Second, there is the bending and twisting 
energy which results from the plate buckling out of its plane. 

A basic assumption is that the plate deflection, W(x,y), is 


infinitesimally small, and that the mid-plane strains may be re- 


presented by 


EE 
e= SZ + + (19) 


er dy tax T Dx oy 


Terms of higher order in W than these are neglected, since W 
is infinitesimal; those terms which are retained are essential for 
the calculation of bending energy. 

In the above expressions Әм, etc., represent the stretching 

ow Y 
of the plate in its own plane. The terms AE RON 5 etc., 
2 \ 3x 
represent the difference between a length of flat plate, and the 


corresponding straight line length when the plate is bent (difference 


between arc and chord). 
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Stresses in the x and y directionsare composed of 

(a) Membrane stresses of finite magnitude, which come 
from the Kirsch Solution and are denoted by Uy, Uy and 

Try below; and (b) Infinitesimal bending and twisting 
stresses which do not appear explicitly below, but which 
are involved in the bending and twisting energy. 


The normal forces per unit length E can be determined from 


N = Th , etc. The membrane energy in the plate is 


X 
| 


Substituting the strains (eqns. 19), апа letting 


UMEM = Ui + U, (21) 
where 
Ag au dv Ou , dv 
„=> o BE + o, BF + (+ BE] de dy (22) 
е т: ЭИ. ЭМ а д 
[ЕЕ +0 84] * xy 9x Bg |^* FY (23) 
an expression for U is found in terms of the stresses, strains, 


MEM 


and the deflection. 
It can be seen that the membrane energy is composed of two 


terms. Omne, U,, relates the work done by the stresses, Ux , etc, 


I 
through a distance сч , etc. The other term, U»: relates the 
oW 
work done by the stresses, Ox, etc., through a distance ف‎ 


etc., as the plate bends out of plane. 

At this point, the assumption is made that the Kirsch 
solution is still valid, even if the plate has buckled infinitesimally. 
This means that U, is independent of the lateral deflection, W, and can 


I 


be treated as a constant. This assumption is justified by considering 


23: 


that the problem is one of incipient buckling, not postbuckling. 
This assumption considerably simplifies the problem, because it 
uncouples the membrane stress distribution from the buckled 
deflection. 


is 


An expression for the bending and twisting energy, 02, 


obtained by considering the moments in the plate and the respective 
rotations at a point as functions of plate deflection. Timoshenko (9) 


shows this expression to be; 


=> b: d - (2-0 DN. Zw 
ES ox? aya 


Adopting the notation М, = = ; Wg7 Sw, etc., and transforming 


e 
zz | dx dy (24) 





the previous energy equations to cylindrical coordinates, with the 
proper limits of integration for the problem of a holé in an infinite 


plate, equations 22, 23, and 24 become, respectively 


U. = CONSTANT (25) 
i h Om Ws 
u, = Bf T [Wye + Je We , 2 Tre We We LE (26) 
4 gg ۳ 
р eT oo 2 
U, = 2f | Wert e ws — (e—9) ws, | + Wee (27) 
O l 





2 i 
— Wee| | rdrdé 
Ү { 
Or Vo and Te are obtained from the Kirsch solution (eqns. 1,2, and 3). 


W and its derivatives are obtained from the assumed series (eq 18). 
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5. SOLUTION OF THE PROBLEM 

With the expression for total energy (UL =U + U, +0.) 
formulated, the next step is to substitute the Kirsch solution, 
and the assumed series for deflection and its derivatives. Then 
this expression must be minimized with respect to each of the 


unknown coefficients, AM or B The substitutions and inte- 


NM" 
grations are shown in Appendices B, C, and D, It is sufficient 

for our purpose here to say that what arises from the integration 
of the expressions for U, and U, are two quadruple series, with 
indices I and K summed over the range of N, and J and L summed 


over the range of M. The quadruple series stems from the multi- 


* 
plication of two double sums. Equations 26 and 27 now become 


Us 5А) ) ) Озь," Аз: "Ак. (28) 
r'y K L 

3 D.) » 2.2. Р.А Axe (2105 
I1 K L 


3 
where the flexural rigidity of the plate, D, equals a 


с 
I 


От ткт, апа PrJKL are different for each value of I, J, K, or L. 


S is the stress applied to the plate at infinity and is factored 
from the Kirsch solution after those expressions are substituted 
into the energy integrals. Ауу апа Аку) are the unknown coefficients 


in the series expression for deflection (eq. 18). 


Following the Rayleigh-Ritz procedure, ш is differentiated 


with respect to each of the unknown coefficients, so that the 


* pate š ; 2 
The unknown coefficients for the cosine series, , only will be 
used henceforth, although these could be the sine series coefficients, 


Bm’ equally as well. 


2 


following set of simultaneous equations result | 


eet DEG Quee tZ 2a) Anat Eu +2 Pa] 


o- Bm = 2 a Qa kt Tt Z M ii Eas, +2 Par A 


O 
" 


(30) 


The last term in each equation arises from the fact that 
BAS (Ar; Are) S 

Аз а practical matter, the equations arrived at in this 
investigation are too complicated to permit considering the 
use of a doubly infinite set of terms in the double series. 
Accordingly, it is henceforth presumed that truncated series 
are used and that there is but a finite number of unknown Am ® 

Simplification in notation is made at this point if the 
quadruple subscripts are replaced by a double subscript arrange- 
ment, 


Let B equal the number of terms taken in M, and & equal 


the number of terms taken in N. Define 
Т'=(т-!)В 5 
J=(K-)B+L 


Using these equations a unique pair of subscripts can be changed 


to a unique single subscript, so that Ату Ата, and Өт ткт, = Өттү! 
етс 
For example, if 4 terms were taken in M, and we wished to 


convert Q to some doubly subscripted Qr: then I = 7, J = á, 


7492 Zu 


К = 9, and L = 2; and 


= (7-1) 4 + 4, a 
= 28 = 34 


so that Q7 492 becomes 28) (34) Similarly 
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A, E A, Gum am Qu 


An ES A, Qu zx Qa 
Аз => A4 О шг > О 
etc. etc. 

Also, let 


J = (TBE = х. В 
where J is the largest value ¿Of I' and J' respectively. 


Substituting 


and 
[on = Чр and Tr's' — Pots if n = ae 
Rr = 2° Огу and pes! = 2° ех if т! > mi 


equation 30 simplifies to 








ОСЕ = 6S (1-0? 
BA ^ 23 | D Ris: + т.) 
qu 
OU = 0= E 
aA, J | EH Ros” + Tay (31) 
т" 
US рези 
A +, 
a 
When the substitution 
_ 6S(I-v) (32) 
À EH 
is made, and the above set of equations expressed in matrix 
notation, the following matrix eigenvalue problem results 
(RA+I)A=o (33) 


where 


29; 


Ra Raa Ro3 
De |" а | (34) 
| Rs, za ze Кр 
IEEE Tr 
Ta Ta ۰ ° چو‎ 
n = т » (35) 
Ta Tar c ° "334 
and 
A, 
Az 
AT EE (36) 
As 


The matrix T is positive definite, since it arose from the 
bending energy integral. The R matrix arose from the membrane 
energy integral, which contains terms from the Kirsch solution. 
This matrix is not necessarily positive definite. Both matrices 
are symmetric. These properties are important to the method 
utilized in solving the matrix eigenvalue problem (see Appendix D). 

The question of whether differentiating ps and setting these 
derivatives equal to zero, produces a minimum or a maximum, has not 
been discussed. Either situation can occur in the problem as it 
has been developed. Some of the eigenvalues obtained in a solution 
of equation 33 represent maxima which are not relevant to the plate 


buckling problem. The criterion for rejecting such answers is the 


set of equations = 
o Ut 

> О 
OA? 





30; 





=> O 
JS 
; (37) 
U > 
on: 


If the set of equations 37 is satisfied, then UL is minimized, 
When the set of equations 31 is again differentiated with 


respect to the appropriate coefficient, the set of equations 37 








become 
O* Ue AR 
т ut Ty > O 
дА? 
2. 
дА% BÉ zs (38) 
e e Í 
= Вет + Те = о 
on 
where Кот! апа Tra are the diagonal elements of the matrices 


В апа Т, respectively. All of these terms are positive numbers, 


since they arise from the "squares" in the energy integrals (i.e., 
no cross product terms). 


By rearranging equation 38 so that it becomes 


Тиз ' = 
AR IU S (39) 
33) 


It is seen that for some negative values of А, and all positive 
values, the solution of equation 33 represents a minimization of 
energy. 

We are interested in only positive values, so that for any 
answer of interest, a minimization is assured. It should be noted 


that the negative values of À that are admissible stem from 


DTE 


buckling modes in the plate that correspond to an applied compressive 
loading. These modes arise from the reversal of sign in the Kirsch 
solution and the associated reversed stress distribution around the 
hole. Modes of this type should be expected from a theoretical 
standpoint. 

When a nontrivial solution has been found for equation 33, 
giving an eigenvalue À and a corresponding eigenvector À, then the 
buckling stress may be found by use of equation 32 and the buckled 


mode shape from equation 18. 
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6. RESULTS AND CONCLUSIONS 

The results of this thesis may be presented compactly in the 
form of a formula for incipient buckling, derived from equation 37, 
viz. 

2 
See = «E )40( 

where ра is the value of 5, the tension at infinity, which causes 
incipient buckling, t is plate thickness, a is the hole radius and 
E is Young's modulus of elasticity. Values of the buckling coef- 
ficient, K, are given in table l. 

Values of K were obtained from a digital computer solution (see 


Appendix F), for the particular case 


il 


30 x 10% psi (Young's Modulus) 
QUE NEC (Poisson's ratio) 


h 


.001 hole radii 

Equation 40 permits prediction of а for other values of t/a 
and E, However, the effect of different values of » cannot be studied 
easily since Y appears in the bending energy integral and in the 
boundary conditions. All results given are for)= .3, but the 
computer program can be modified by changing card listed as 0011 
on page 84. 

Four distinct mode shapes are apparent, as determined by 
different buckling coefficients appearing for each. A tabulation 


of K is shown in Table 1 below. 
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Table 1 


Values of the Buckling Coefficient, K 


MODE K 
1 755695 
2 2.830 
3 Liny 
4 12.3 


Since the energy method involves the use of a series which 
cannot be made infinite in the present method, the actual buckling 
coefficients obtained will always be higher than the true value. 
This is because the deflection of the plate is only approximated, 
and the energy never quite minimized, A solution with less plate 
energy could be found, if more terms were taken in the series for 
deflection. 

Different values of the buckling coefficient, K, actually 
calculated by use of the digital computer program are shown in 
Table 2, and are plotted against the number of terms taken in 
both the N and M parameters of the deflection series, equation 18. 
These plots are shown in Figures 4a through 4e. The projected 
values listed in Table 1 are taken graphically from these figures 
as the asymptotes which the curves for buckling coefficient appear 
to approach, when the number of terms in N and M increase. These 
projected values are considered to be the best evaluation of the 
buckling coefficient. 

Figure 4f is taken from the experimental data of Danis (14). 


All of his points for a finite plate are plotted as K versus d/b, 
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Table 2 


Values of the Buckling Coefficient, K, for 


Numbers of Terms in the Series 


9.986897 
921333 
8.9990 
8,9927 
8.9730 
8.9347 
8.8897 


9,1250 
8.2133 
8.4310 
8,4253 
8.4023 
8.3637 
8.3253 


16.288 
15.0959 
15.910 
152307 
15,882 


14.761 
14.463 
14.415 
14.412 
14.392 


Terms in M 


8.6493 
8.0673 
7.9743 
7.9683 
2.99597 
7. 029 
7.99 


6 


8.4420 
2.9967 
7.8060 
Dr 997 
7.7847 


Terms in M 


63587 / 
8.0630 
7.9860 
2.9773 
29577 
109251 
7.8940 


6 


8.2107 
99 
129197 
7.9100 
7.8907 


Terms in М 


L52592 
527] 
12.679 
12.516 
12.493 
12.488 
12.470 


14.545 
12.616 
12.046 
112897 
115572 


Terms in M 


16.418 
13.714 
122929 
12,719 
12.682 
172077 
12-662 


19.709 
13.318 
12.589 
127388 
12250 


954 


200 
„173 
‚8763 
„7857 
mu 7 


м Jj N سم من‎ 


3033 
913 
>50 
29050 


- Y Y 00 


14.371 
I SOIT 


1107095 


15.660 
1955259 


120957 


Various 


8.4157 
7.8748 
7.7840 
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8.5026 
7.9906 
7.9107 


14.348 
12.487 
11.0929 


15.647 
297 
125599 


8.4153 
7.8746 


8.5026 
749906 


14.347 
12.486 


15.646 
L52252 


the ratio of hole diameter to plate width. The buckling coefficients 
in Table 1, for modes 1 and 2, are plotted on the y axis (for zero d/b). 

The values for modes 1 and 2 are in fair agreement with experi- 
mental data. No correlation between the theoretical prediction and 
the experimental observation is attempted for modes 3 and 4, 

The buckled deflection of the plate for modes 1 through 4 is 
shown in Figures 5 through 8 respectively. It is seen in these 
figures that the four possible conditions of symmetry are satisfied 


by the four modes. This result is given in Table 3. 


Table 3 


Resultant Symmetry for the Various Mode Shapes 


MODE SYMMETRY WITH RESPECT SYMMETRY WITH RESPECT 
TO X AXIS TO Y AXIS 
1 £ (y)we s&(-y) (х= (зл) 
2 OS TEY) TEES) 
3 t(y) 23E) CC) tee (GX) 
- (у) = ==) i in 


In looking at Figure 4&, it is apparent that mode 2 is 
converging to a value of buckling coefficient that is somewhat 
higher than that for mode 1. This can be qualitatively justified 
by considering the symmetry of buckled shape, and the equivalent 
one dimensional problem. 

Consider two beams, buckled by some unspecified loading, into 
shapes as shown in Figure 9. 


EL A {e co 


a b 


Figure 9 


DIE 


It is known from beam theory that the shape of Figure 9b 
has a higher critical buckling load than the shape of Figure 9а. 
This is exactly the relationship between modes 1 and 2 of the 
buckled plate, as can be seen in Figures 5 and 6, respectively. 
There is a small area near the intersection of the hole 
edge and the y axis, which assumes cross sectional shapes Similar 
to that in Figure 9a for mode 1, and in Figure 9b for mode 2. 
It is only in these small areas connecting the "top" of the plate 
to the "bottom", that this sort of bending occurs, so that the 
disparity between the buckling coefficients of the first and 
second plate modes is much smaller than that for the beam modes. 
This same argument applies to the difference between modes 
3 and 4, and to the gross differences between modes 3 and 4, and 
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- Plot of the Buckling Coefficient 
Mfor mode 1. 
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for mode 2. 
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APPENDIX A 
SATISFACTION OF THE BOUNDARY CONDITIONS 
The free edge boundary conditions for the hole are stated 


in equations 15 and 16 as 





W We , (2-0) 3-0 = 
Were $ Wee — We + ES”) Wee es Wee = O (A-1) 
r=) 
and 
W ee A-2 
[wet о (үу, + =] = о (А-2) 
r=) 
The plate deflection, W, is taken as 
оо оо 
_ A... cos M6 
w(ne)=) у. r "fi> Хым + Yun | NM oy (A-3) 
fe Bam SIN MO 


where the coefficients Акм апа ВАМ are varied to minimize the plate 


energy, and the coefficients X Mm and Yum are chosen so that the 
boundary conditions expressed by equations A-1 and A-2 are satisfied. 
The derivatives of W are tabulated in Table 4. When these 


derivatives are substituted into equations A-1 and A-2, and evaluated 


at r= 1, the following equationsresult 








Anm cos МӨ 
- 3 [s + 5, Хум + 53 Үмм RU o 
N=| M=0 Bum SIN MO 
and 
2 е А мм COS MO 
SS eee ftn a 
Bum SIN МӨ 
N*| M=O 
where 
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S = N(N+I) - N3 - MV (A-6) 


5, = (№+1)(№+2) – (+1) 0 – М) 


(А-7) 

S, = (N+2)(N+3)-(N+2)D - MD (A-8) 

£,= —N(N+I)(N+2)+N(N+1) +N + (2--0)MÍN + (3-0) mM" (A-9) 
t,=-(N+1)(N+2)N+3)+(N+1)(0+2) + (N+1) + (2-D) M”(N+1) + (3-0) M (С=О; 
t¿=- (N+2) (N+3)(N+4) + (N+2)(N+3) +(N+2)+ (2-2) M'(N*-2) «(3 -) M* (A-11) 


Equations A-4 and A-5 must hold for ail values of 0, which can 


only be true if 


S, + S; Хам + Sa Yan A (A-12) 
+ ES + P. =0 (A- 13) 


for any value of N and M. Thus, a value of X Mm and YNM can be 


found for each value of N and M, such that the boundary conditions 


are satisfied. 


By application of Cramer's Rule: 

















=S, 53 
Xn = ima з = Ta + 5з t, 
M 
S, S; S203 7 S405 
to ta 
and 
S, -$, | 
vA = tz ps = “ Sat: + Sı to 
j Sa S3 52 Ез ж Sato 
t2 ts 








о 


The above results are incorporated into SUBROUTINE BOUND (see 
APPENDIX F) which generates values of XM and Y M for given values 


of N and M. 
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APPENDIX B 


EVALUATION OF THE ENERGY INTEGRALS 


The total energy in the plate, from equations 25, 26, and 27, is: 


where 


U, = constant 


bf flew? + Sows . 200 We 4 (B-1) 
О аз [Wr + + E dr 49 
О l 


2 


D 2'T r oo 
= >| E Wer Ма. 21 (B-2) 
© l 


201-0) We ur Leej- пе Mol” r dr 49 


= 
| 


Note that these integrals contain second degree terms in the 
derivatives of the deflection, W. Since W is represented by a 
series, care must be taken to insure that the series are proper ly 
multiplied. 

The product of two double sums, both summed over the same range 
for the subscripts N and M, is a quadruple series, with dummy subscripts 
I and K summed over the range of N, and K and L, summed over the range 
of M. 


For example; 


co oo e OO oo сю oo 
) ) tier” )cos M6| = У) (17) соз гө ) (lr) cos LO 
N=] M=0 J=; J=0 k=1 L=0 


оо 
3 (1+2) (1+7) соѕ Je cos L© 


1=0 


Me 
1 
Ns 


н 
CA 
" 
© 
х 
M 


D Mer ra [e E Cos JO cos LO 


NU 
[x] 
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It is seen that the cross product terms above do not combine 
: : f ~I 

directly into one term. One cross product term involves г ¿the other 
| БЕ | . 
involves r . For this reason, cross product terms in the energy 
integrals B-1 and B-2 must be identified and separated. 

As an aid in keeping the cross product terms separated, the 
notation iS adopted that starred (*) functions are functions of I 
and J, and that non-starred functions are functions of K and L, in 


place of N and M, respectively. Equation B-1 may be written 


p | [ww + теша. 
о 


і 











(B-3) 
+ +T (WW, + W.w2)| r dr d8 
and with the substitution of the Kirsch solution becomes 
v (1-2 = АМА. +01 3, - Z|WEW. cos 26 (B-4) 
s - m om cos 20 
Equation B-2 is now 
D 2T. oo 
ولا‎ = a J (Usa + Use + Use| rdr dé (B-5) 
о i 
where 
= 4 
U a = (Wry + Xe + Woe | UM + Wr + Wee (B-6) 
x 
Us, 7 73 Ju [ore SE] | + “| (8-7) 
* 
Us auc) - a Wes - we (B-8) 
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A listing of the expansion of the products of the derivatives 
necessary for the evaluation of U, and U, are shown in Tables 5 


and 6, respectively. Further presentation of the energy integrals is 
in the notation delineated in those tables. 


Equations B-4 and B-5 now become 


U,= УУУ р-з, COS EOTE TEE (B-9) 


- (Е, +36, cos 20 -|Е,,+Е,,- Ea ЗЕ КЕСЕ 2Egp] Sin 20] de 


љ= 25у) у | ЧГ +в Gat Ga Gg (B-10) 
О 


+ бу. + Gg +(1-0) [Gon + Gop+ бз + Gor] 


 2(-3)|6; * Gs + баа + Sa, | de 


When the integration. in 9 is done, equations B-9 and B-10 reduce to 
oo ec oo Da 
1 Arskı E (B-11) 
T=! J=0 k=) 50 Bry’ Bri 


and 


Us 3325 Baki a (B- 12) 


T=1 $20 k-11=0 p e e. 
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Tables 5 and 6 


Tabulation of the Integrals of the Radius Dependent Terms 


NOTATION: 


G ...ور‎ :. eC. are functions ot IA 


1? da Eos 
x means X Mm as used previously, only as a function of K & L. 
x” means X AM as used previously, only as a function of I and J. 
y means Yam as used previously, only as a function of K and L. 


y^ means Y as used previously, only as a function of I and J. 


NM 


THE INTEGRALS ARE TABULATED ON THE FOLLOWING PAGES. THE FIRST 
NINE TERMS IN EACH COLUMN ARE SUMMED TO FORM THE WHOLE EXPRESSION. 
THE TERMS a(J,L) or b(J,L) SHOULD BE CHOSEN FOR THE UNKNOWN 


COEFFICIENTS Акм or Bm RESPECTIVELY. 
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Tabulation of the Compacted Form of the Terms in the 
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Table 5 


Membrane Energy Integral 
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Table 5 (continued) 


59 


M. i IJ "KL: 


tal | 
b(J,D) M 


ub Ld “Ay Agp 
"6b b(3,L)- тта Bic, | 
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b(J,L) ‘Arg ‘Axr 


Table 5 (continued) 
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where 


a(J,L) соз ЈӨ соз LO 


b(J,L) = sin JO sin LO 


IK 
TISE42) 


M I(K*l)x 
VR 
, I(K-42 

x: B 
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1+K+3 
‚ (1+1) (К+1)хх“ 
rm 
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1+K+ 
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^ `° 


(I+2)K 
T+K 


n (тг) (Кэз )жу” 
+K+ 

' (142) (eee yy 
I+K+ 


where 
a(J,L) = cos J@ cos LO 
b(J,L) = sin J@ sin LO 


Table 5 (continued) 
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where where 
a(J,L) = cos J@ cos LO a(J,L) = sin J0 sin LO 
DTD) =з sin LO b(J,L) = cos JO cos LO 
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Table 5 (continued) 
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where where 
a(J,L) = sin JO sin LO a(J,L) = sin Je sin Le 
b(J,L) = соз ЧӨ соз ГӨ b(J,L) = cos JO cos LO 
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Table 5 (continued) 
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Table 5 (continued) 
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Table 6 


Tabulation of the Compacted Form of 


the Bending Energy Integral 
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where where 
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where where 
a(J,L) = cos J@ cos Le a(J,L) = cos Je cos LO 
b(J,L) = sin J® sin LO b(J,L) = sin Je sin LOG 


Te 


where 
a(J,L) 
b(J,L) 


II 


Table 


IJKL 
(I+K+2) 


тло 
T +K+3 
IJ(K+2)L 
Lo 2 


(ї+1) кіх“ 
+К+З 


(1+1)J(K+1)Lxx" 
Кау T 


(1+1) J(K+2)Lyx" 
T+K+ 


(I-2)JKLy 


T+K+ 
(I*2)J(K41)Lxy^ 
T+K+ 
(I*2)J(K*2)Lyy^ 
T+K+ 


sin Je sin LO 


cos JO cos LO 


6 (continued) 
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where 
a(J,L) 
b(J,L) = 


13 


JK 


~ (I+K+2) 


J(K+1)x 
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J(K+2) 
IFR 


tà 


ТЮЗ 
J(K+1)xx 
+K+ 
J(K+2) x* 
“тл 


JK 
I+K 


J(K+l) s 
"TORRES 
Tar 

1+K 


sin JO sin LO 


cos JO cos LOG 


Table 6 (continued) 
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where where 
a(J,L) = sin JO sin LO a(J,L) = sin JO sin LO 
b(J,L) = cos JO cos LO b(J,L) = cos JO cos LO 
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APPENDIX С 


INTEGRATION OF SINE-COSINE TRIPLE PRODUCTS 


When evaluating the energy integrals for the plate, it is 


necessary to evaluate the following two integrals 


2 
L= Í [cos AO cos BO cos ce|de (Eso) 
O 


and 


217 
ie =f AO sın BƏ cos co) de (Ger) 
о 


where 4, B, and C are integers. 
The above integrals can best be evaluated by expanding the 


integrands using the identities 


SIN x cos u = > [їн (х+у)+ SIN (-9)] (C-3) 
COS "X SIN = El (x+y) -SIN (x-y)| (C-4) 
Cos X cosy= [cos (x+y) + cos(x-y)| (C-5) 
SIN X SIN y = 4 |cos(x-y)- cos(x+y)| (C-6) 


The object is to reduce the integrals to the form 
e "IT 21 
[sm DO de ov | De de (C-7) 
0 0 
where D is an integer sum or difference of À, B, and C. Note that the 
sine integral in C-7 is always zero and that the cosine integral is 


nonzero only if D is zero. 


Expanding equation C-2, using identity C-3, we get 


27 
I Jl SIN Ao|[si(8+c)O + siN(B-C)0| d6 (C-8) 
0 


yo 


Continuing the process with identity C-6 


aT 
Is = + [с°$(А-8-с)ө- сов (А+®+с)ө +cos(h-8+c)0-cos(A+8-c)0 [do (C-9) 
0 


I. is now expressed as the sum of four integrals 


І;= 1+1, +13+ 14 (C-10) 
where 
T/2 if (A-B-C) = 0 
L = 
O otherwise 
-"I/2 if (A4B4C) 2 0 
L, = 
O otherwise 
11/2 if (A-B4C) - 0 
I, = 
O otherwise 
-T/2 if (A4B-C) = 0 
I, = 


O otherwise 


А similar expansion for I. (equation C-1) can be made. The 
results are the same as above, except that the signs of I, and L, 


are positive. 


76. 


APPENDIX D 


SOLUTION OF THE EIGENVALUE PROBLEM (Rx + T) A = 0 
The canonical form of the matrix eigenvalue problem is 


(1A-B)X=0 0-1) 
where I is the identity matrix, X an unknown eigenvector, А 1з 
an unknown scalar parameter (eigenvalue), and B a known coefficient 
matrix. 

Equation D-1 bas a solution if determinant of 

(IN-B)=0 (D-2) 
If this determinant is expanded, a polynomial equation in N results. 
However, if the determinant is of order greater than three or four, 
it is generally inefficient to obtain this polynomial equation and 
then solve it. 

There are alternative procedures for solution of the eigenvalue 
problem. A method which gives extremely accurate eigenvectors and 
eigenvalues for symmetric matrices is one named for Jacobi, which 
is presented by Fox (1). Brown and White (13) give a comparison of 
the accuracy of the Jacobi method with that of other methods available. 
The necessity for extremely accurate eigenvectors is shown later. 


Consider the eigenvalue problem (equation 33) 


[(RAFT)A =0 (D-3) 


where both R and T are symmetric matrices and T is positive definite. 


Dividing equation D-3 by A , and making the substitution 


| 
Б (D-4) 
gives 


(ly -R)A=0 (D-5) 


I: 


Any symmetric matrix with real elements can be reduced to a 


product of the eigenvector and eigenvalue matrices, so that 


li 
NO 


* 
Ba (D-6) 
where D is a diagonal matrix of the eigenvalues (4,) ot pondo 
is a matrix of normalized eigenvectors CHF obtained from a 


solution of the problem 


ien 


1-41 


\ai= б (D-7) 


т = : Е . . 
Note that Q = 0, since the matrix T is symmetric and has orthogonal 
eigenvectors. 


Equation D-5 may now be written 


(QDQy-R)A=0 (D-8) 
Premultiplying by Qu. we get 
(29y-QR)A=0 (D-9) 


5 = Т И 
Factoring a Q out of A = ТА = Q Q À, and letting 


Y=Q'A (D-10) 
the problem reduces to 
(ps-g'&G)Y «c о 


Since T is positive definite, D contains all positive terms. 


By again premultiplying, this time by D 


na 


NY 


, factoring D ^, and making 


< 
= 


the substitution 


zeap*y (D-12) 


equation D-11 reduces to 


т 


| =! 
(z4-g^g'ReE *)z -o (0-13) 
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and becomes 
(Iy-8)2=0 (D-14) 


which is of canonical form. 


L ES 
Note that D^ and D ^ are easily obtained by taking the square 


roots, and the square root reciprocals, respectively, of the diagonal 
terms. Also note that B is symmetric. The problem is now well 
suited for solution by use of the Jacobi method. It is seen that 
inaccuracies in obtaining eigenvectors Q would cause compounded 


errors in the final solution for A and À 2 


The problem is returned to original form with the substitutions 


Az- (D-15) 


M 
3 


and 


(D- 16) 
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APPENDIX E 


SUGGESTIONS FOR EXTENSION OF THE PROBLEM 


There are two obvious directions for extension of this problem 
of a thin plate with a circular hole. The buckling problem for 
a strip of finite width and the postbuckling problem for both finite 
and infinite plates remain to be treated analytically. 

The postbuckling problem is the more difficult of the two 
because of the coupling between the stress distribution and the 
buckled deflection. Johns (2) outlines the postbuckling situation 
for plates and suggests several methods of solution. Solution of 
the postbuckling problem would be important because of the revised 
stress concentration factors that would be obtained; cf. Houghton 
and Rothwell (16). 

Extension of this problem to a finite strip, of width 2b, would 
require the following: 

a. Using the Howland solution (17) instead of that due to Kirsch, 

to obtain the stress distribution in the plate. 

b. Satisfying additional boundary conditions for two "free" 

straight edges. 

Integration of the energy terms for the finite strip is between 


b 


the limits r =4 and r = Ena > so that the integration with 
respect to O becomes more complicated than for the infinite plate. 
This integration is over the interval 0 to 27. Wylie (10) 
presents a method for the evaluation of real trigonometric integrals 
over this interval. The method involves substitution of appropriate 


complex variables for the sine and cosine terms, and then evaluation 


of the residues inside a unit circle. 
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There are two possible methods for satisfaction of the 
additional boundary conditions of the finite strip. Introduction 
of additional unknowns, similar to Xm and YM in equation A-3, 
into the assumed series for deflection could be made. These 
unknowns could be determined by simultaneous solution of all the 
boundary condition equations. Another method, the Lagrangian 
multiplier method, introduces dummy variables so that additional 
restraining equations, those expressing the boundary conditions, 
can be solved simultaneously with the energy expressions. 
Budiansky and Hu (12) present a clear example of the use of this 
method in minimizing the plate energy, while satisfying plate 


boundary conditions. 
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APPENDIX Е 
COMPUTER PROGRAM LISTINGS 


The following pages contain a listing of the computer programs 
used in the solution of the problem. All programs are written in 
FORTRAN IV for use on an IBM System 360 computer. 

The control card listing on the first page is effective as of 
September 1967. Users of this program should check current procedures. 

The maximum terms to be taken in N and M, and the plate thickness 

"read into the computer by data card (see comment cards regarding 
input requirements). Young's modulus, E, and Poisson's ratio are 
established by cards serial 0010 and 0011 in the main program. 

Printed output from the computer is a column matrix EIVU which 
is an array of all of the values of S cr in psi for the parameters 
of plate thickness, Poisson's ratio, Young's modulus, terms in N, 
and terms in М, Matrix EIVU is followed by a map of the buckled 
plate deflection for the minimum value of tensile critical buckling 


stress. 
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